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ABSTRACT: 

An  analysis  of  a  3- inch  gun  launched  finned  motor  case  is  performed  in  an 
attempt  to  determine  the  cause  of  failure.   Two  specific  features  of  the  problem 
are  identified  as  the  probable  reasons  for  the  buckling  failure  of  the  motor  case. 
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INTRODUCTION 

The  Naval  Weapons  Center,  at  China  Lake,  California,  is  presently 
in  the  process  of  developing  a  series  of  gun  launched  guided  projectiles. 
Recent  experimental  firings  of  a  3 -inch  projectile  showed  the  structural 
design  of  the  motor  case  to  be  inadequate  since  the  case  buckled  severely 
as  the  projectile  passed  through  the  gun  barrel.  The  authors  were  asked 
to  perform  an  analysis  of  the  motor  case  to  determine  the  cause  of 
failure.  In  addition,  the  authors  were  asked  to  recommend  a  structural 
design  of  a  5  inch  projectile  that  minimizes  the  motor  case  weight  and 
to  propose  a  design  technique  which  can  be  used  by  engineers  as  a 
design  tool  in  future  gun  launched  guided  projectile  projects.   This 
report  concentrates  on  the  analysis  of  the  3-inch  motor  case  and  the 
design  technique. 

The  following  four  pages  show  l)  a  drawing  of  the  projectile; 
2)  and  3)  high  speed  photographs  of  two  firings;  and  k)     a  picture  of 
a  buckled  projectile  in  flight. 
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DESCRIPTION  OF  THE  THREE -INCH  PROJECTILE 

The  casing  is  Ul30  steel  alloy  heat  treated  to  200,000  psi.  From 
MTL-RDBK-5A,  Feb  8,  1966: 

Elastic  modulus  in  tension  and  compression,  E  =  29  x  10  psi 

Modulus  of  rigidity,  G  =  11  x  10  psi 
Therefore 

Poisson's  ratio,  v  =  .32 

Density,  cu  =  O.283  #/in3 

0.002  Compressive  yield  stress,  a   =  198,000  psi 


Total  projectile  weight: 
Warhead  with  filler 
Fuse 

Motor  case 
Propellant 
Forward  bulkhead 
Nozzle  plate 
Fins 

Fin  Release 
Obturator 
Nozzle  plug 

Length ,  L  =  9 . 5  in 
Radius,  r  =  1.5  in 
Thickness,  h  =  0.08  in 


19.80  lbs 


forward  =  11.1*3  lbs 


motor  =  5.58  lbs 


aft  =2.79  lbs 


DESCRIPTION  OF  THE  EXTERNAL  LOAD 

The  external  load  on  the  projectile  consists  of  two  parts;  l)  the 
breech  pressure  at  the  base  of  the  projectile,  and  2)  a  lateral  pressure 
on  the  sides  of  the  case  due  to  the  breech  pressure  seeping  beyond  the 
obturator.  A  typical  plot  of  the  predicted  barrel  pressure  as  a  function 
of  time  is  shown  in  Fig.  1.  According  to  NWC  (Data  sheet  for  5  inch 
projectile,  L.  R.,  10-20-71)  the  base  pressure  on  the  projectile  is 
0.9^-2^  x  breech  pressure.   Thus,  a   ,  the  predicted  maximum  rigid  body 

ttlSIX 

acceleration  of  the  projectile,  is 

maximum  P 

.  force  _  (26.^  ksi) (0.9^2*0 (tt) (1-5  in.)  _  fl 
amax  "  total  ~  19.8  lbs/g  "  °'yUU  g 

mass 

where  g  is  the  acceleration  due  to  gravity. 

The  lateral  pressure  in  the  motor  case  occurs  when  the  obturator 

cracks .  Photographs  of  the  firings  show  smoke  coming  from  the  gun 

barrel  before  the  projectile  emerges,  thus  establishing  the  existence 

of  this  pressure.  No  estimates  of  the  magnitude  of  this  pressure  are 

available.  Future  NWC  obturator  designs  will  attempt  to  eliminate  this 

pressure. 

DETERMINATION  OF  THE  INTERNAL  AXIAL  FORCE 

In  order  to  obtain  an  estimate  of  the  transient  internal  axial  force 
in  the  motor  case,  the  projectile  subjected  to  the  transient  base  pressure 
shown  in  Figure  1  is  idealized  as  a  one  dimensional  steel  structure  with 

Tit 

a  force  P  sin  ' "  ,  0  £  t  £  .012,  at  the  aft  end.   The  forward  end  of 
o     . 012 ' 

the  projectile  is  taken  as  either  free  or  fixed,  these  two  conditions 


being  the  extreme  possibilities.   The  belief  was  that  this  analysis 
would  provide  a  bounded  estimate  of  the  transient  axial  force.   The 
details  of  the  analysis  and  the  results  are  presented  in  Appendix  A. 
The  essential  result  of  the  analysis  is  that  the  rise  time  of  the 
pressure  is  very  slow  in  relation  to  the  wave  speed  of  the  case,  and 
hence  the  internal  axial  force  at  any  location  along  the  case  at  any 
time  is  equal  to  the  mass  forward  of  that  location  multiplied  by  the 
rigid  body  acceleration  of  the  projectile  at  that  time.   Thus,  there 
is  neither  an  amplification  nor  a  reduction  in  the  internal  axial 
force  due  to  dynamic  effects. 

Since  the  axial  force  in  the  case  is  dependent  upon  the  distribution 
of  the  mass  of  the  projectile,  the  question  of  the  load  path  of  the 
motor  propellant  becomes  important.  The  propellant  body  force  could 
either  be  transferred  to  a  longitudinal  shear  load  on  the  case  or  a 
compression  load  on  the  forward  bulkhead  of  the  aft  end,  depending 
to  some  extent  upon  the  bonding  condition  between  the  propellant  and 
the  case.  Since  there  was  considerable  uncertainty  as  to  which  path 
was  the  correct  one,  a  two  dimensional  stress  analysis  was  performed 
on  the  propellant  for  both  the  condition  of  bonding  and  of  no  bonding. 
The  Rohm  and  Hass  computer  program  AMG032A  for  the  stress  analysis  of 
linear  elastic,  isotropic,  axisymmetrically  loaded,  axisymmetric 
propellant  gain  with  (or  without)  an  orthotropic  elastic  motor  case 
was  used  (Ref.  l).   The  details  and  results  are  presented  in  Appendix 
B.  The  essential  conclusion  drawn  from  the  results  for  the  bonded 
propellant  is  that  the  maximum  shear  stress  on  the  case  (and  in  the 
bond)  is  approximately  300  psi  and  approximately  80$  of  the  total 


VS.  "T  i  *>£ 


propellant  body  force  is  carried  by  longitudinal  shear  loading  on  the 
case  and  20$  by  compression  on  the  bulkhead  at  the  aft  end  of  the 
case.   Of  course,  the  unbonded  propellant  body  force  is  carried 
entirely  by  a  compression  on  the  bulkhead. 

Taking  the  results  of  these  analyses  into  consideration,  the  forces 
at  the  forward  and  aft  ends  of  the  motor  case  are 

forward  =  (^f^)  (8>9°°  g)  =  102  kips 

Paft  =  102  kips  +  (1-99  +  g8  X  3^9  lbs)(8,900  g) 

=  1^5  kips 
for  the  bonded  propellant  (80$  carried  by  the  case)  and 

Paft  =  102  kips  +  (^r)(8>900  g)  =  120  kips 
for  the  unbonded  propellant.  The  maximum  stress  due  to  these  loads  is 


CTmax  =  (lU5  kiPs)/C(2Tr)(l.5  in.)(0.08  in.)] 
=  192  ksi 

which  is  just  below  a 

cy 


DETERMINATION  OF  THE  DYNAMIC  RESPONSE  OF  THE  CASE  TO  THE  TRANSIENT  AXIAL 

FORCE 

Because  the  axial  force  is  on  the  case  for  only  a  short  period  of 
time,  the  possibility  exists  that  the  case  may  not  suffer  large  deforma- 
tions for  loads  that  are  larger  than  the  critical  load  for  static 
buckling,  i.e.,  the  loading  may  be  impulsive.  This  possibility  is 
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investigated  in  Appendix  C,  where  the  critical  static  buckling  mode  of 

nt 
the  case  is  subjected  to  the  force  P  sin  n  p  and  the  elastic 

transient  response  is  computed.   Two  estimates  of  the  natural  frequency 

of  the  critical  buckling  mode  are  used,  one  without  the  propellant 

mass  attached  and  one  with  the  propellant  mass  attached.  The  essential 

result  from  this  analysis  is  that  if  the  maximum  value  of  the  applied 

force,  P  ,  is  equal  to  the  critical  load  for  static  buckling  the 

maximum  magnitude  of  the  buckled  deflection  is  approximately  100-200 

times  larger  than  the  initial  imperfection  in  the  shape  of  the 

critical  buckling  mode.  For  the  manufacturing  tolerances  in  this 

program  this  is  a  large  deflection  and  represents  a  severe  buckling 

condition.  Thus,  the  loading  is  not  impulsive  and  the  maximum  axial 

force  must  be  less  than  the  static  buckling  load  to  prevent  severe 

buckling . 

DETERMINATION  OF  THE  STATIC  BUCKLING  LOADS 
OF  THE  CASE 

The  results  of  the  preceding  analyses  indicate  that  the  motor 
case  can  be  designed  using  standard  procedures  for  estimating  the 
static  buckling  strength  of  shells.  Consequently,  the  following  static 
loading  conditions  were  considered: 

(la)  Uniform  axial  compression,  only,  P  =  Pporwa  d  =  paf+ 

(lb)  Linearly  varying  axial  compression  only,  P  =  Pforwara  +  tK3  in)(x)q£ 


A  value  of  P  equal  to  90%  of  the  critical  load  for  static  buckling 
causes  an  amplification  of  approximately  10  which  may  also  represent 
severe  buckling.  This  depends  upon  the  magnitude  of  the  imperfection 
that  exists  in  the  buckling  mode  shape  prior  to  firing. 
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(2)  Uniform  external  radial  pressure  only,  q 

(3)  Combined  uniform  axial  compression  and  uniform  radial 

pressure,  P  =  P  _     ,  =  P  „.  ,  q 
^       '       forward    aft'  Hr 

where  q  and  q  are  the  radial  and  longitudinal  uniform  external 

pressure  components  respectively,  and  x  is  the  longitudinal  coordinate 

(x=0  at  the  forward  end). 

The  elastic  buckling  load  for  loading  condition  (la)  was  determined 

using  the  formula  presented  in  Ref .  2,  page  528,  and  the  computer  program  for 

the  geometrically  nonlinear  analysis  of  arbitrarily  loaded  shells 

of  revolution  developed  by  Ball  (Ref.  3).  The  results  are  shown 

in  Fig.  2,  where  a   is  the  critical  axial  buckling  stress.  The 

D   '        cr 

maximum  predicted  axial  stress  in  the  case  due  to  the  firing  is  192 

ksi  and  is  also  indicated  in  Fig.  2.  Since  a       is  much  larger  than 

cr  to 

a  ,  the  buckling  is  inelastic  and  the  predicted  elastic  buckling 
loads  are  invalid.  The  inelastic  buckling  load  for  loading  condition 
(la)  was  estimated  using  the  charts  from  Ref.  2,  page  704 ,  and  is  also  given 
in  Fig.  2.   Note  that  this  stress  is  less  than  the  maximum  predicted 
axial  stress  in  the  shell.  As  a  consequence  of  this,  the  case  will 
buckle  inelastically  due  to  the  axial  load  alone. 

The  elastic  buckling  load  for  loading  condition  (lb)  was  determined 
using  the  computer  program  of  Ref.  3>  and  is  shown  in  Fig.  2.  For 
this  condition  all  the  propellant  was  assumed  to  be  carried  by  a 
uniformly  distributed  shear  stress  on  the  motor  case.  The  ratio  of 


-* 
The  result  shown  in  Fig.  2  is  a  conservative  estimate  since  NWC  experimental 

results  for  the  static  buckling  of  a  O.065  in.  thick  cylinder  indicated 

a  higher  buckling  stress. 
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axial  stress  at  the  forward  end  to  the  stress  at  the  aft  end  was 

■   *j  o  =  .67^}  the  same  ratio  as  the  total  mass  ahead  of  the 
forward  end  to  the  total  mass  ahead  of  the  aft  end  of  the  case.  The 
critical  stress  shown  in  Fig.  2  is  the  stress  at  the  aft  end. 

Buckling  due  to  loading  condition  (2)  is  elastic.  The  results 
from  Ref.  2,  page  5U0,  and  the  computer  program  (Ref .  3)  are  presented 
in  Fig.  2,  where  q   is  the  critical  external  radial  pressure. 

For  elastic  buckling  due  to  loading  condition  (3),  Ref.  2,  page 
5^U,  recommends 

*+-SL-l 

P     q 
cr   ^-cr 

The  result  from  the  computer  program  is  given  in  Fig.  2.  Note  that 
this  result  agrees  with  the  recommendation  of  Ref.  2.  No  results 
for  the  inelastic  effects  for  loading  condition  (3)  are  readily- 
available  .  A  computer  program  that  can  treat  inelastic  buckling  of 
cylinders  under  general  loading  must  be  used  to  obtain  the  inelastic 
buckling  load  under  that  loading  condition.  Hence,  the  stability  of 
the  shell  along  the  dashed  line  is  unknown  at  the  present  time. 
Also  shown  in  Fig.  2  is  the  maximum  axial  stress  possible  i.e. 
200,000  ksi.  Hence,  the  cross-hatching  in  Fig.  2  indicates  the 
limiting  values  of  the  allowable  loading.  From  this  it  can  be  seen 
that  if  the  lateral  pressure,  due  to  leakage  beyond  the  obturator,  is 
larger  than  2800  psi  the  case  will  buckle. 

The  critical  elastic  buckling  mode  predicted  by  the  computer 

TT  X 

program  for  all  four  loading  conditions  was  essentially  sin  -=-  cos  28 

jj 

where  9  is  the  circumferential  coordinate. 
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The  details  of  the  analyses  using  the  formulas  and  charts  in 
Ref.  2  are  given  in  Appendix  D.  The  computer  results  have  been  given 
to  NWC  personnel,  Code  U573. 

An  examination  of  Fig.  2  reveals  that  there  are  significant 
differences  in  the  predicted  elastic  buckling  load  obtained  from  Ref.  2 
and  from  the  computer  program  for  loading  conditions  (la)  and  (2). 
Apparently,  there  is  need  for  further  study  here.  The  shell  considered 
here  is  not  exactly  a  thin  shell,  and  the  formulas  in  Ref.  2  may  not 
be  valid  for  this  shell. 

DESIGN  TECHNIQUE 

The  following  plan  is  recommended  as  a  general  design  technique: 

1.  Determine  the  loads  P_     , ,  P  _, ,  a  ,  and  q  from  the  estimated 

forward   aft'  rr'     ^s 

weights  and  acceleration  data  and  a  propellant  stress  analysis. 

2.  Select  the  thickness  of  the  case  such  that  the  maximum  stress  is 
less  than  198  ksi. 

3.  Check  to  see  if  P  _.  and  q  need  to  be  changed. 
J  aft     us  to 

k.     Determine  if  the  shell  is  stable  (inelasticially)  under  the  loads 

P.o     j  5  P  ia »  q  s  q  •  Note  that  q  should  include  both  the  effect 
forward'  aft'  Hr'  ^s  ^r 

of  the  slumping  propellant  and  the  effect  of  the  external  pressure. 
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CONCLUSIONS 

Several  analyses  have  been  performed  on  the  3-inch  motor  case  to 
determine  the  cause  of  the  buckling  failure.  The  results  of  these 
analyses  are: 

1.  The  dynamic  effects  are  insignificant  and  the  case  responds  to  the 
breech  pressure  as  if  the  pressure  were  applied  in  a  static 
sense,  i.e.  the  maximum  internal  axial  load  in  the  case  can 

be  computed  on  the  basis  of  the  mass  distribution  and  the 
maximum  rigid  body  acceleration  of  the  projectile. 

2.  The  loading  is  not  impulsive,  i.e.  it  is  on  the  case  for  a 
sufficient  length  of  time  such  that  the  shell  will  collapse  if 
the  applied  load  equals  the  static  buckling  load  of  the  case. 

3.  Approximately  80$  of  the  propellant  inertial  load  is  carried 
by  shear  along  the  case  when  the  propellant  is  bonded. 

k.     The  case  will  buckle  inelastically  due  to  the  axial  load  caused 
by  the  launch. 

5.  If  the  obturator  cracks,  10-20$  of  the  breech  pressure  acting 
as  a  lateral  pressure  is  sufficient  by  itself  to  cause  elastic 
buckling. 

6.  The  predicted  buckling  mode  for  axial  load  and  for  lateral 
pressure  is  the  same  and  is  essentially  identical  to  the 
observed  buckled  shape  of  the  launched  projectile. 

7.  The  static  buckling  loads  computed  using  the  formulas  in 
Ref .  2  are  considerably  lower  than  those  predicted  by  the 
computer  program  of  Ref.  3« 
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APPENDIX  A 
Transient  Axial  Response 

Statement  of  the  Problem 

In  this  section,  the  case  is  idealized  as  a  one  dimensional, 

nt 
one  material  bar  with  an  input  force  function  P(t)  =  P  sin  —  , 

^  o     a 

where  a   is  the  duration  of  the  force  and  P  the  maximum  intensity 
at  t  =  tt/2.  The  time  length  of  the  input  force  was  taken  as  .012  sec. 
The  boundary  condition  at  the  aft  end  of  the  bar  (x=0)  was  taken 
as  free,  that  is  the  axial  force  in  the  bar  is  equal  to  the  forcing 
function.  Two  boundary  conditions  were  considered  at  the  forward 
end  of  the  bar  (x=£); 

i)   a  free  end  with  axial  force  equal  to  zero,  and 
ii)  a  fixed  end  with  displacement  equal  to  zero. 
The  actual  boundary  condition  at  the  forward  end  is  somewhere  between 
these  two  idealized  conditions.  Hence  the  actual  solution  lies 
between  the  bounds  obtained  from  the  free -free  and  free-fixed  problems 

The  physical  parameters  taken  in  this  problem,  were  the  length 
of  the  bar  I  =   10  in . ,  and  the  Young ' s  modulus  for  a  steel  bar 
E  =  30  x  10  lbs /in.  . 


P(t) 


aft 


forward 


A-l 


Governing  equation  for  axial  response 


B2u  r2     S2u     J^-       for  moving  bar 

2        2  ~  ]  W 

St        Sx      I   0     for  stationary  bar 


where  u  =  displacement  along  x  axis,  x  axis  attached  to  the  bar 
t  =  time 

c  =  speed  of  wave  propagation  =  y—  (2) 

E  =  Young's  modulus  of  elasticity 
p  =  mass  density 
I   =  length 
A  =  area 

The  initial  conditions  are 

u   (x,0)  =  0  (3) 

ut  (x,0)  =  0  (10 

where  u,  denotes  velocity   ^r-  . 

t  0"C 

The  boundary  conditions  are 
i)  for  free-free  boundaries 

u  (0,t)  =  G(t)  (5) 

ux  (J,t)  =  o  (6) 

ii)  for  free-fixed  boundaries 

ux  (0,t)  =  G(t)  (7) 

u  (i,t)  =  0  (8) 


In  equations  (5)  through  (9) 


u   denotes  -^— 
x  Sx 


P(t) 
where  A  is  the  area  of  the  bar. 


G(t)  =  S|i  (9) 


A- 2 


I.  Axial  response  with  free-free  boundary  conditions  with 

P(t)  =  P  sin  2*.  (10) 

o      a  v 

Then 

G(t)  =-°-sin^=  P  sin  ^  (11) 

We  seek  solution  of  equation  (l),  with  initial  conditions  (3)  and  (h) , 
and  boundary  conditions  (5)  and  (6),  with  forcing  function  given  by 
(10). 

Solution  by  superposition.  Consider  solution  u(x,t)  in  the  form 

u(x,t)  =  v(x,t)  +  U(x,t)  (12) 

Then  (l)  becomes 

vtt +  «tt  - c2  *j  +  <W  ■  -  fii  -  -2 1 sin  ^       (13) 


Letting 


f(x,t)   .  -  Utt  +  c2  U^  -  |M  W 


equation  (13)  becomes 


\t  -   °2  v«  =  f(x't}  (15) 


Boundary  conditions  (5)  and  (6)  become 

v  (0,t)  =  G(t)  -  U  (0,t)  (16) 


x 


v   (£,t)  =  -U   (X,t)  (17) 
and  initial  conditions  (3)  and  (k)   become 

v  (x,0)  =  -U  (x,0)  (18) 

vt(x,0)  =  -Ut(x,0)  (19) 
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U(x,t)  is  chosen  so  that  equations  (l6)  and  (17)  become  homogeneous 

boundary  conditions  in  v.  Here  we  take 

2 
U(x,t)  ?  (x  -  |j)  G(t)  (20) 

With  equation  (20)  we  find 

i)  boundary  conditions  (16)  and  (17)  become 

v  (0,t)  =  0  (21) 

v  (l,t)  =  0  (22) 

ii)  initial  conditions  (l8)  and  (19)  become, 

2 
v(x,0)  =  -  U(x,0)  =  (x  -  |-)  G(0) 

But  G(0)  =  0,  hence  v(x,0)  =  0 

2 
vt(x,0)=  -  Ut(x,0)  =  (x  -  §j)  Gt(0)  =  Y(x) 

-      2 
Since  G  (0)  =  P^  Y(x)  =  &  (x  -  |j) 

v(x,0)  =  0  (23) 

vt(x,0)  =  Y(x)  (2k) 

where 

Y(x)  =  ^  (X  "  fl}  <25) 

The  solution  to  the  boundary  value  problem  in  v(x,t)  defined  by  equations 

(15),  (21)  -  (2k)   is  sought.  A  solution  in  the  form  of  a  Fourier 

series  is  assumed,  say 

00 
v(x,t)  =  £  vn  (t)  cos  2J™  (26) 

n=0= 
which  satisfies  boundary  conditions  (21)  and  (22).  Functions  f(x,t) 

and  Y(x)  are  also  expanded  into  their  Fourier  series, 
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where 


f (x,t)  =  )  f  (t)  cos  SJJ* 

n=0 


fo(t)  =  \    J  f(|,t)  d? 

0 


(27) 


(28) 


fn(t)  =|  J  f(5,t)  cos  2Jp  d§    n  > 

0 


(29) 


and 


T(x)  =   £     Tn  cos  ££ 


x 


n=0 


where 


'o-J    I    T<5>d? 


0 


fn  -J     *  (?)  cos  ^d?      n  >0 


0 


(30) 


(31) 


(32) 


Substituting  expressions  (26),  (27)  and  (30)  into  equations  (15)  ? 
(23),  and  (2k)   yields 

£«*¥'{  ^^  +  C'  (  T^  '»<*>  -  *„<*>}  ■  °       (33) 


n=0 


dt 


I  cos  E«  Vn(o)  .  0 


n=0 


£  4  v  (0)     -. 


m 


(35) 


n=0 
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Equations  (33)  -  (35)  are  satisfied  if  for  each  n,  (n=0,l,2, . . . ) 
d2  v  (t)         .2 

— 72-  +  <=  ("t)    *»<*>  "  '*<*>  (36) 

at 
vn(o)  =  0  (37) 

dv  (0) 

l£—  -  '„  <38) 


Equations  (36)  -  (38)  may  be  numerically  solved  on  the  computer  (using 

Runga-Kutta  or  some  other  scheme)  for  v  (t),  for  n=0,l,2,...,  where 

the  tilda  denotes  the  numerical  approximation  solution.  Hence  for  a 

convergent  series, 

00  00 

v(x,t)  =Y    vn(t)cos  ^f  ^Ivn(t)cos  2J"U*(x,t)  (39) 

n=0  n=0 

Once  v(x,t)  is  obtained,  we  return  to  equation  (12)  to  find 

u(x,t)  =  v(x,t)  +  U(x,t)  °»  v(x,t)  +  U(x,t)  =  u(x,t)  (U0) 

Hence , 

N 

u(x,t)  -  £  ^n(t)  C0S  *1T  +  U(x't}  {kl) 


n=0 
The  force  is  given  by 


N 


F  =  EA  fc  =  {-  7    2LH   v  (t)  sin  ^  +  g  (x't}j  EA  (1*2) 

dx        l    Z-      j£       nv  £         dx  J 

n=l 

Note  that  F  does  not  depend  upon  the  n=0,  or  rigid  body  term. 
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Numerical  solution  of 


<*2  v  (t)  .    ,      x2 

nv  2   /iitt\  /.n 
^ +  c      (—  )  v   (t) 

dt2           u  ■■  n 


f  (t) 

n 


n  =  1,2,.. 


(36) 


v  (0)   =  0 


(37) 


dv  (0) 
n 

dt 


=  ¥_ 


n 


(38) 


Recall  eq.    (29)  f or  n  >  0 


f,(t)-|  J  *ct.t)  eo.  24s  «- §/[(-?  ♦£)<** 


0 


d§ 


2  To       J         1       _  nng         1     -     .     nng  2g  n  tt  g 

=  -r      G,  .     ■{-   -. rn-  COS   r*  "  "/ T   §   sin  a       +  /      \o    cos   « 


2JL 


2  .     n-nt    L         £ 

! TTT      Sin  y*  +  r*- 

nn\3  H 


I  J 


tt  G  + 

X  plA.  /    n-n 


P    N     I      .     ring 
—  sin  — -* 


Jo 


sin 


nrrll 
4     J 


Therefore 


or 


where 


f  (t)  = 

nv 

21 

2     2     tt' 
n     tt 

21 

2     2 
a    n 

f     =  K 
n        n 

.        TTt 

sin  — 

K     =  - 

n 

2J&  P 
2     2 

=.        .        TTt 

P  sin  

a 


(MO 
(U6) 


A-7 


V 

Let  p  =  =r-  ,  then  eqs .  (36)  to  (38)  become, 
n 


&2   P  (t)    0  ,  '  v2  , 


"^2-  "  +  C  W     Pn  =  Sin  T  (W 

Pn(o)  =  0  (k8) 

dp       Y 

-w  (0)  ■  r  <■*> 

n 
Recall  equations  (32)  and  (25), 

0 

y     =  2      P     ttP    (-       i\   CQS  n_n_i  _  2_rrP  /-£\  ,      , 

Yn       £     J        a     U       2i}  cos       X  ai    l2     2J  ^0) 


n     rr 


Y 

-£-  =  +  SL  =  +  3.82  x  10"3  (51) 

Kn 

for  a  =   .012  sec. 

Hence  we  solve  for  Pn(t)  from  the  following  equations 

§ +  =2  (xf  \  ■  ^  *r  <*> 

dt 

Pn(0)  =  0  (53) 

dp  (0) 

_JL —  =  *  =  3.82  x  10_:5    for  a  =  .012  (5*0 

CL"C       TT 

After  obtaining  p  (t)  from  the  computer  we  find 

v  (t)  =  K  P  =  -  -~-^  X  P  P  (t)  (55) 

n      n  n     2  2      n  v  ' 

n  a 
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Finally  then  the  force  per  unit  P  is 
N 


nn  Vn^   .  nnx  L  1  3U(x,t)  \  raJ1 


p  *     p 

n-l 


Therefore 

N 

F     V  2  tt   Q  /.\   .   n  ttx   /..   Xn   .   irt  /_^-v 

-  =  I    —    P  (t)  Bin  —  +  (1  -  j)  sin  —  (56) 

°   n=l  a  n 

The  latter  term  is  the  internal  force  due  to  the  rigid  body  acceleration. 
The  terms  under  the  summation  sign  introduce  the  transient  dynamics 
due  to  the  elasticity  of  the  bar. 

II.  Axial  response  with  free-fixed  boundary  conditions  with 

P(t)  =  Pq  sin  ^  (10) 

A  solution  of  equation  (l),  with  initial  conditions  (3)  and  {k) ,  and 
boundary  conditions  (7)  and  (8),  with  forcing  function  (10)  is  sought. 
The  method  proceeds  as  before,  except  here  we  take 

U(x,t)  =  (x-1)  G(t)  (57) 

Boundary  conditions  on  v(x,t)  become 

vx(0,t)  =  0  (58) 

v   (A,t)  =  0  (59) 

The  governing  equation  in  v(x,t)  is 

vtt  -  c2  Vxx  =  f (x,t)  (60) 
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where 


f(x,t)  .  -  Utt  +  c2  U^  (61) 

=  -  (x-X)  Gtt 


The  initial  conditions  are 


where 


v(x,0)  =  0  (62) 

vt(x,0)  =  (x-i)  Gt(0)  =  Y(x)  (63) 


Y(x)  =  -^  (x-X)  (6U) 

a 

Assume  the  Fourier  expansion 

00 

v(x,t)  =  £  vn(t)  cos  ^JL*  (65) 

n=l,3,... 

which  satisfies  boundary  conditions  (58)  and  (59)'  As  before  f(x,t) 

and  Y(x)  are  expanded  into  Fourier  series, 

00 
f(x,t)  =    £      fn(t)  cos  2JL*  (66) 

n=l,3,... 
2X 
f0(t)     =^1    I    f(S,t)  d§  (67) 

0 
24 

f  (t)     =i     f    f(5,t)  cos  2JLld?        n>0  (68) 


V   '         X 


2X 


y(x)  =  Y    \  cos  ^r  (69) 


n=l,3,... 

24 


Yn  =  T    J       Y(?)   COS  ^    d?  n  >  °  (70) 


0 

2X 


\  -  h  !    *<?>  cos  ^  d? 
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Substituting  eqs.  (65)  -  (70)  into  eq.  (60),  (62)  and  (63)  yields 

no  "it 

00 

I   cos  S^E  vn(0)  =  0  (72) 

n=l,3, .. . 

£  r  dv  (0)      <i 

I  -^  {^T--\}  =  °  (73) 

n=l,3,... 


Again  we  seek  the  solution  to  the  equations 
2 

=2(fi)Vf„  <*) 


dt2 

vn(0)  =  0  (75) 

dv  (0) 

~^~  -  \  (76) 

for  n  =  1,3,5,...   Numerical  solution  for  a  finite  number  of  n  =  1,3,5,. 

N  yields  the  approximate  solution,  v  (t),  hence 

N 
v(x,t)  =  £   vn(t)  cos  £|T£  (77) 

n=l,3,... 

Then 

N 

u(x,t)  -  £   vn(t)  cos  2J^£  +  (X.X)  G(t)  (?8) 

n=l,3,-.. 

Finally 

N 

|(X,t)  -IS  V*)  si"  TFT  +  G<*>  (79) 

n=l,3,. . . 
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Numerical  solution  of 


a2v 


dt 


n  .  2 
2"  +  c 


iitt  \ 
2  A  / 


v  =  f 
n    n 


v  (0)  =  0 
n 


dv  (0) 
n 

dt 


=  Y_ 


n 


for  n  =  1,3,5, 


Calculation  of  f  (t)  and  Y  for  n  >  0 
nv       n 

24 


21 


fn^  -  T  /  f(?'t}  C°S  T^  d?  =  T  I  (^  Gtt  COS  ^  d?   <*» 


G^"  fn(t) 
Gtt  n 


0 


2     2 

n  rr 


for  n  =  2,1+, 6,.  . 
for  n  =1,3,5,. • 


(81) 


^-JJ«0~^«-i  J  ^(M)eoS^id§ 


0 


(82) 


Tn   .012 


0    for  n  even 

~8^  * 

for  n  odd 


2  2 
^  n  tt 


(83) 


Hence  we  seek  the  solution  of 


d2yn(t)    2  Ait/      -8£P 
5—  +  c   5T)  v  =  5- 


dt 


v  (0)  =  0 

n 


V2X/   n 


.     TT  t 

2    2     .012 
a   n 


dv 


dt 


*(0)  = 


■81 


2     2 

n  rr 


ILL 

Qi 


for  n=l,3,5, 


(8U) 

(85) 
(86) 


Let   K 


n 


= s-  ,  and  define  3  =  v  /K  , 

2   2'  n    n  n  ' 

a       n 


(87) 
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then  the  above  equations  become, 


d2P 


dt2 


"  +  "   ©  *n  "  sin  "4  (88) 


MO)  =0  n  =  1,3,5,...  (89) 

E  (o)  =  -  =  +  3-82  x  10'3  for  a   =  .012  sec  (90) 


n 

dB 


dt   V  '     TT 


After  6  is  determined,  we  have 

n  ' 


v  (t)  =  K    B     =       ~8  V  2     Mt)         n=l,3,5...  (91) 

n  n     n        (.012r  n         n 


and 

N 

F  V  Utt      „   /.x      .     n  ttx         .       rrt        ,~ns 

-  =    I  —$-  pn(t)  sin  -3j-  +  sin  -y-      (92) 

F°  n=l,3,...   *  n 

for  the  solution  of  a  free-fixed  bar.  The  latter  term  represents  the 
static  response  to  the  end  force.  The  terms  under  the  summation 
sign  introduce  the  dynamic  effects. 

III.  Results 

The  time  history  F/P  at  x/4  =  0,0.25,0.5,0.75  and  1  are  given 
in  Figs.  A1-A5  for  the  free-free  bar  with  a  =  .012  sec.  Figs.  A6-A10 
show  the  corresponding  values  of  F/P  for  the  free-fixed  bar.  Figure 
All  shows  the  history  of  F/P  at  x/Jl   =  0.5  for  a  =  0.0012  and  Fig.  A12 
shows  the  history  at  x/JL   =  0.5  for  a  =  0.012  and  F  =  3  x  10  psi  for 
the  free-free  bar.  These  figures  are  from  a  plotter.  The  solution 
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was  plotted  after  each  new  term  in  the  series  was  added.  The  solution 
with  the  contribution  from  U  was  plotted  last.  Hence  the  small 
wiggles  indicate  the  dynamic  contributions.  Figure  A13  shows  F/P 
as  a  function  of  x/X  for  several  values  of  t  for  a  =  0.012  sec.  for 
both  bars. 
IV.  Summary  of  Results 

The  solutions  of  the  boundary  value  problem  show  that  for  a  forcing 
function  of  time  duration  .012  seconds,  the  force  appears  as  a  standing 
wave.  That  is,  the  time  length  of  the  force  is  so  much  larger  than 
the  time  of  wave  propagation  that  the  dynamic  effect  is  negligible. 
Decreasing  the  time  of  the  force  pulse  to  .0012  seconds  showed  a 
noticeable  increase  in  the  dynamic  effect.  However  even  at  this 
shorter  time  pulse,  the  internal  forces  are  but  slightly  larger  than 
those  for  a   =  .012  seconds.  A  further  decrease  of  <y  to  .00012  seconds 
show  that  the  dynamic  forces  are  the  same  order  as  the  static  result 
and  there  is  an  amplification  of  the  internal  force  in  the  bar. 

With  the  time  pulse  at  a  =  .012  seconds,  Young's  modulus  of  the  bar 

6  6 

was  reduced  from  30  x  10  psi  to  3  x  10  psi  to  determine  the  effect 

of  "softening"  the  bar  material  due  to  inelastic  effects.  It  can  be 

observed  in  Fig.  A12  that  although  a  noticeable  increase  in  the  dynamic 

effect  occurs  the  internal  force  is  not  significantly  altered. 

The  conclusion  that  the  dynamic  effects  are  insignificant  can  be 

supported  by  an  examination  of  equation  (52).  For  each  value  of  n, 

equation  (52)  is  equivalent  to  a  single  degree  of  freedom  system  with 

a  harmonic  forcing  function  of  frequency  rr/a  rad/sec  or  260  rad/sec 
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for  a  =  0.012  sec.  The  natural  frequency  of  the  single  degree  of 
freedom  system  is  cnrr/i.  For  i  =   10  and  c  =  2.02  x  10  in/sec,  the 
lowest  natural  frequency  is  63,500  rad/sec.  Hence,  the  ratio  of  the 
forcing  frequency  to  the  lowest  natural  frequency  is  0.00^1.  The 
magnification  factor  for  this  frequency  ratio  is  essentially  one, 
i.e.  the  case  response  in  a  static  sense. 
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APPENDIX  B 
Two-Dimensional  Stress  Analysis  of  the  Motor  Propellant 

The  results  for  the  propellant-case  interface  stresses  from  the 
Rohm  and  Hass  computer  program  AMG032A  (Ref .  l)  are  given  in  Fig.  B-l 
for  E  =  5  x  10  psi  and  a  body  force  of  ^80  g,  which  is  equivalent  to 
an  acceleration  of  8000  g  of  the  propellant  grain  whose  density  is 
0.06  lb/in  .   Also  shown  in  Fig.  B-l  is  the  location  and  the  magnitude 
of  the  maximum  principal  tensile  stress.  In  this  analysis  the  motor 
case  was  assumed  to  be  rigid.  The  nodes  of  the  finite  elements  used 
to  model  the  propellant  are  denoted  in  Fig.  B-l  by  I,  J,  where  I 
defines  the  radial  location  and  J  the  axial  location  as  shown  in  Fig. 
B-l.  The  stresses  predicted  by  the  program  are  at  the  center  of  each 
element.  Hence,  the  stresses  shown  in  Fig.  B-l  are  actually  located 
0.125  in.  from  the  propellant-case  interface.  A  plot  of  the  shear 
stress  from  the  inner  face  of  the  propellant  to  the  outer  face  is 
shown  in  Fig.  B-2  for  J  =  lU.  Note  that  the  outer  face  propellant 
shear  stress  is  approximately  10$  larger  than  the  stress  given  in 
Fig.  B-l  at  J  =  Ik, 

The  total  vertical  force  provided  by  the  longitudinal  shear  stress 
along  the  case  and  the  axial  compression  stress  on  the  bulkhead  at 
the  aft  end  is 


The  computer  program  output  consists  of  the  coordinates  and  displacements 
of  each  nodal  point  and  the  stresses  and  strains  for  each  continuum 
element . 
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-Coo     -4fog>     '2co 


T~j[     Z     3    f  .f       "^   ..-*»  ...- 


^ift.    B-l        Prop«IUrt-f-Stressc5    a,+    Jnfe- 
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Reaction  =  (°'25  ±n'^2^   +  2?3  +  •-.  58  +  19  psi)(l.l)(Tr)(3  in.) 

+  (0.25  in.)[(930  psi)(l.25  in.)  +  (787  psi)(l.75  in.) 

+  (832  psi)(2.25  in.)  +  (731  psi)(2.75  in.)psi](TT) 

=  22.9  kips  (shear)  +  5.0  kips  (compression)  =  27.9  kips 

when  the  shear  stresses  shown  in  Fig.  B-l  are  increased  by  10$  to 
extrapolate  to  the  interface  stresses .  The  total  force  applied  to 
the  propellant  is 

Papplied  =  (0-°6  lb/in.3)(l.52  -  0.52)(in.2)(9.5  in.  )(tt) (8,000  g) 
=28.6  kips 

Thus,  the  stress  analysis  appears  to  provide  total  equilibrium.  The 

fact  that  P    ,  .   is  slightly  less  than  P   .,  .  ,  is  due  to  the  fact 
reaction  applied 

that  the  compression  stresses  at  the  propellant -bulkhead  interface 
are  somewhat  larger  than  the  values  shown  in  Fig.  B-l,  which  are 

at  the  center  of  the  bottom  elements  0.125  in.  above  the  bulkhead. 

22  Q 
From  these  results  on ,':?  =  82$  of  the  total  propellant  inertial 

load  is  carried  by  longitudinal  shear  stress  on  the  motor  case. 

k 

The  same  problem  was  analyzed  with  E  =  5  x  10  psi.  The  results 

were  essentially  identical. 

The  stresses  for  the  unbonded  propellant  with  E  =  5  x  10  psi  are 
shown  in  Fig.  B-3.  The  total  vertical  force  is  carried  entirely  by 
compression  on  the  bulkhead  at  the  aft  end  and  is  equal  to 

Reaction  =  (°*25  in)[(^-l8  psi)(l.25  in.)  +  (4,^90  psi)(l.75  in.) 
+  (4,516  psi)(2.25  in.)  +  (4,526  psi)(2.75  in.)](Tr) 
=28.2  kips 
The  computer  output  for  these  three  runs  have  been  given  to  NWC 
personnel,  Code  4573* 
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APPENDIX  C 
Transient  Lateral  Motion 

Statement  of  the  Problem 

The  object  of  this  investigation  was  to  determine  the  transient 
lateral  behavior  of  an  axially  loaded  bar  with  initial  deformation 
w  (x).   The  bar  is  idealized  as  simply  supported  at  both  ends,  and 

is  subjected  to  an  axial  compression  P(t)  =  P  sin  where  a  =  .012 

seconds.  The  behavior  of  the  bar  was  considered  for  various  ratios 
of  P  j  the  maximum  intensity  of  thrust,  to  P  ,  the  Euler  buckling 
load  for  a  simply  supported  column. 


C-l 


,vtA 


± 


P(t) 


/yf^ 


\y}  w 


.___..  J» 


B.C.'s  simply  supported 


w  -  initial  imperfection,  assume  w 


w  -  bending  deformation,  assume  w 


-  y 


n=l 


-I 


1      .     n  rrx 
w     sm  — t— 
n  I 


x,   t.\      .      n  ttx 
fn(t)   sm  — r 


(1) 


(2) 


n=l 


o  J _. 


v 
h* 


/ra.0      .JU     =  inertia  force 

i  ^     d 
w  =  w     +  w 


y '    M.    =  0:      M  -    (M  +  dM)   +  S  (m  ^~)   +   (V  +  dV)   dx  +  P  dw 
u      A  dt 


=  0 


_|M+V  +  Pdw  =  0 
dx  dx 


.2 

d  w 


)   V  =  0:      V  -    (V  +  dV)   -  m  2*  dx  =  0 

u  at* 


dV  B2w 


(3) 
GO 
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Recall  the  moment  curvature  relation, 


M(x)  =  -  EI  ¥C 


xx 


(5) 


Using  (3)  and  (5),  equation  (U)  becomes 

ft  5  w    -p  9  w  -  ™  S  w 
-EI  jj-  -  P  — £-  =  m  — 2" 

hx  Bx      St 


(6) 


With  (l),  (2)  and  (7),  equation  (6)  becomes 


"  I    EI  <!f>k   *n<*>  si*  ^T  +  ^    I    <T>2  ("n  +  fn<*»  si«  ^ 


n=l 


n=l 


=  m 


I  *»W  Sin  —i 


(8) 


n=l 


which  may  be  written 


V  [ei  (Sf)    fn(t)  -  p(t)  £f)2  (w*  t  fn(t))  +  m  r<t)~ 

n=l 

n  =  1,2, . . . 
Equation  (9)  is  satisfied  if  each  of  the  coefficients  of  sin 


n  ttx   _ 
sin  — t—  =  0 


(9) 


n  ttx 


k  2 

EI  (2?)  f  (t)  -  P(t)  Fp)      (wi  +  f  (t))  +  m  f'(t)  =  0         (10) 


I  '       n 


n    n 


n 


with  I.C.'s 


w(x,0)  =  w  (x) 
dw 


5t 


(x,0)  =  v(x) 


(11) 
(12) 


C-3 


Using 


-p 

= 

2 

n 

2 

TT 

EI 

e 

/ 

% 

= 

2     2 

n     tt 

i2 

i 

EI 

m 


(13) 


(HO 


equation  (10)  becomes 


?.(t)  ♦  I  a  -  ^i)  f„(t) .  £M  .Si 


n 


n 


pe 


n  n 


n  =  1,2,... 


with  initial  conditions 


f  (0)  =  w  X 
n      n 

f  (0)  =  0 
nv  ' 


where 


w 


n 


i  /  >.   .   n  ttx  , 
w  (xj  sm  — t—  dx 


(15) 


(16) 
(17) 

(18) 


For  an  impulse  load  in  the  form, 


p<*>  ■  po sin  -m 


and  introducing  the  amplification  factor 


f(t) 


n 

the  previous  system  of  equations  become, 

P  .  - 


••/,\         2  f.          o.  irt  "1               o.        rrt          rc 

gn(t)  +  con  [1  -  -  sm  -^J  gn  =rsm—     con 

e  e 

gn(0)    -  1 


gn(o)  =  o 


(19) 


(20) 


(21) 
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This  equation  is  also  valid  for  a  cylindrical  shell  subjected  to 

an  axial  load  or  an  external  pressure  where  to  is  the  natural  frequency 

in  the  n   mode  and  P  is  the  buckling  load  in  the  n   mode. 

e 

The  previous  system  was  solved  numerically  on  the  computer  using 

a  Runga-Kutta  integration  method  to  yield  the  results  shown  in  Figs. 

C1-C9.  Several  values  of  P  /P  were  considered  for  both  to  =  6,000  rad/sec 

o'  e  n  ' 

and  10,000  rad/sec. 
Summary  of  Results 


Amplification  Factor  |g(t)  | 

THEIX 

p  /p 

o'  e 

to  =  6000  rad/sec 

to  =  10,000  rad/sec 

.75 

3 

k 

.90 

10 

10 

1.0 

100 

180 

1.1 

6oo 

220,000 

1.2 

^50,000 

co* 

1.3 

.* 

00 

J.  TO 

oo  defined  as  greater  than  1  x  10   [maximum  g(t)  possible 
from  computer  integration  algorithm.] 

Estimation  of  the  Natural  Frequency 


K  a  Gh  ..        i?x  h2 
n   aim      K2'  a2 


m  =  mass /unit  area 


From  "Tables  of  Frequencies  and  Modes  of  Free  Vibration  of  Infinitely 
Long  Cylindrical  Shells,"  by  M.  Barron  and  H.  Bleich,  J.  Appl.  Mech, 
Vol.  21,  195^,  PP.  178-18U. 
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where 

a  =  1.5  in . 

G  =  11  x  103  psi 


,-h   #-sec  ..  oQ  .     _  -,,   ,  -5  #-sec 


m  =  7.*+3  x  10   SZEjp  x  .08  in.  =  5-9^  x  10"^  Z-~ 


in  in* 


when  the  mass  and  elasticity  of  the  propellant  is  neglected, 

For  w  =  sin  —7-  sin  26*     (—  =  6) 
*  a 

K  =  0.09178         ~  =  2331 

IT 


and  hence 


10  =  9«6  x  10  rad/sec 


Adding  the  mass  of  the  propellant  gives 


m  m  5.9U  x  10"5  +  I.38  x  10   x  tt  x  1.52/(3tt)  =  1.62  x  10   Z"l£2. 

in 


Hence 


u)  =  5.8  x  10  rad/sec 
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APPENDIX  D 
Computation  of  the  Static  Buckling  Loads 

Axial  Compression  Only  (Ref.  2,  page  528) 


acr  _  n     Eh 
~T  "     c  r 


where  acr  is  the  design  -  allowable  buckling  stress,  i.e.,  nine  out  of 
ten  shells  will  not  buckle  at  a  .  For  elastic  buckling  T]  =  1.  From 
Fig.  3-23-1,  page  530  of  Ref.  2. 

C  =  0.23 
c 


^E-cTolh18-8*  and 


z      ^(l-v2)1/2     (9.5)2(l-;322)1/2   _. 
Z  "     rh      _      (1.5)(0.08)   "  715 


Therefore 


a 
cr 


=  (0.23)  (29xl03  ksi)  (0.08  in.)/(l.5  in.) 


=  356  ksi 


which  is  larger  than  a     .  Hence  the  buckling  is  inelastic,  i.e. 

cy 

T|  <  1. 


This  value  of  r/h  is  beyond  the  end  of  the  curve,  and  hence  the  value 
selected  for  Cc  is  an  estimated  one.  This  indicates  that  perhaps  this 
value  is  invalid  for  this  small  r/h  ratio. 
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Figure  3.62-8,  page  70U  of  Ref.  2,  gives  the  critical  inelastic 

buckling  stress  o"    as  a  function  of  a  /71  for  elastic  buckling 
cr  cr'  '  to 

for  alloy  steel  Ul30,  heat  treated  to  180  ksi.  From  curve  E  , 

a   =  162  ksi 
cr 

This  stress  is  based  upon  a  0.002  compressive  yield  stress  of  16U.5  ksi, 

which  is  considerably  below  the  stresses  of  179  ksi  and  198  ksi  given 

in  MTL-HDBK-5A  and  for  U130  alloy  steel  heat  treated  to  180  ksi  and 

200  ksi  respectively.   Increasing  the  predicted  critical  inelastic 

buckling  stress  to  account  for  the  higher  heat  treatment  leads  to 

acr  ~  (162  ksi)  (i2|)  =  179  ksi 

when  the  MIL-HDBK  stresses  are  used. 

A  static  buckling  test  was  conducted  at  NWC  on  an  axially  loaded 
3  in.  cylinder  of  U130  alloy  steel,  heat  treated  with  h  =  O.O65  in. 
For  this  cylinder, 

a 

-§r  =  290  ksi 

Tl 

and  hence  from  Fig.  3*62-8  of  Ref.  2  the  inelastic  buckling  stress  is 

a       =   159  ksi 
cr 

so  that 

acr  «  (159  ksi)  (i||)  =  176  ksi 

From  the  experimental  results,  the  axial  force  at  buckling  was  121.3 
kips .  Hence 

(rr         \  -  121.3     kjpS  .     lq?r     y-i 

^cAxp  ~   (tt)(3  in.)(0.065  in.)  "  197*5  kS1 
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which  is  considerably  above  the  predicted  value  of  176  ksi  and  is 

essentially  equal  to  a  .   This  very  high  value  is  somewhat  surprising, 

cy 

Nevertheless j  in  light  of  this  the  predicted  inelastic  buckling  load 
of  179  ksi  for  the  0.08  in.  case  must  be  considered  low  (conservative) 
External  Radial  Pressure  Only  (Ref.  2,  page  5*+0) 

%  -  K  ™\      (r)2 

H    Pl2(l-v2)   L 

From  Fig.  3.23-7  of  Ref.  2,  lateral  pressure  only, 


K  =  2k 
P 


for  Z  =  715.   Therefore 


CTcr  =  ,2hj    (tt)2  (29xl03ksi)(0.08  in.)2 
T\  (12)  (1-.322)  (9-5  in.) 

=  U5  ksi 

which  is  below  <j  .  Hence,  this  buckling  is  elastic,  and  the  critical 
lateral  pressure  is 

q   =  !^  =  (U5  ksi)  ^'f^O     =  2.1.0  ksi 
^-cr    r     v         (1.5  in.) 
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